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£NJ ■ Abstract 

Using analytic and numerical techniques Nielsen-Olesen vortices, which in the con- 
text of Ginzburg-Landau theory are known as Abrikosov vortices of type-II supercon- 
ductors, are studied for large Ginzburg-Landau parameter A. We show that their energy 
is equal to (im 2 /2) log A to leading order, where n is the winding number of the vortex, 
(-h ' and find that the limit of the gauge field can be expressed in terms of the modified 

Bessel function K\ . The leading terms of the asymptotic expansion of the solution are 
given, and the different contributions to the energy are analyzed. 

> 

Q\ 1 PACS numbers: 11.15.Kc, 11.15Mc, 11.27+d 
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^ : 1 Introduction 

Of all the localized finite-energy solutions of classical gauge theories, the vortices of the 
Abelian Higgs model in 2 space dimensions, the prototype of a gauge field theory with 
spontaneous symmetry breaking, should be the ones easiest to understand. However, none 
of these solutions is given in terms of known functions. Nielsen and Olesen [T] found 
the time-independent, radially symmetric, localized finite-energy solutions of the Abelian 
Higgs model in 2 space dimensions, the Nielsen-Olesen vortices, by reducing the equations 
of motion to two second-order equations for two radial functions. The mathematically 
rigorous proof that the resulting equations for the two radial functions have solutions with 
the required properties was given by Tyupkin et al. [2J, and Berger and Chen [3J. 

In the context of Ginzburg-Landau theory, which is the time-independent Abelian Higgs 
model without an electric field, the Nielsen-Olesen vortices are known as Abrikosov vortices 
of type-II superconductors [I]. This means that the properties of Nielsen-Olesen vortices 
can, and have been, studied in experiments. The Nielsen-Olesen vortices also provide a 
simple example of cosmic strings [5j , which might explain some of the structures seen in the 
universe today. 
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With the solution not available in terms of known functions, numerical computations 
become all the more important. For the Nielsen-Olesen vortices the numerical work started 
soon after the solutions were found [6] [7J. Asymptotic analysis is another technique often 
applied when the explicit solution is not known. For the Nielsen-Olesen vortices, Berger 
and Chen [3j obtained some asymptotic results for large Ginzburg-Landau parameter. The 
asymptotic analysis of the monopole structure was given by Kirkman and Zachos [8]. 
More recently, the same techniques were used for the Skyrmion [U] and a Skyrme-like 
monopole [10j . In this paper, we perform a similar asymptotic analysis for the Nielsen- 
Olesen vortices. 



2 Radially symmetric solutions 

The Hamiltonian density of the time-independent Abelian Higgs model in 2 space dimen- 
sions is given by 

h = /•;,/•"•' + i(A0(sVr + ^(1 1 2 -i) 2 - (i) 

Here = di<p — lA^ and = d{Aj — djAi (i,j=l,2) are the covariant derivative and 
the field strength, respectively, and the metric is g = diag(l, 1). H in Eq. ([1]) is also the 
Ginzburg-Landau free energy of a superconductor. In this model, the Ginzburg-Landau 
parameter A is equal to 1 at the point between type-I and type-II superconductivity. The 
corresponding Euler-Lagrange equations are 

A£>V-^>(M 2 -i) = o, 

diF** + l - - 0(£>W] = 0. (2) 

The Euler-Lagrange equations have radially symmetric solutions of the form 

4> = f(r)e me , A t = - a ^-e VJ x\ (3) 
where n = ±1, ±2, ... is the winding number. The radial functions satisfy the equations 

a"--a' + / 2 (n-a) = 0, /" + -/'- = ^(/ 2 - 1)/, (4) 

and the boundary conditions for regular vortex solutions to exist are 

/(0) = o(0) = 0, lim f(r) = 1, lim a(r) = n. (5) 

r— >oo r— >oo 

These solutions are the Nielsen-Olesen vortices [1] of the Abelian Higgs model; or for A > 1 
the Abrikosov vortices of type-II superconductors. The existence proof for such solutions 
was given by Tyupkin et al. [2]. The proof is based on the fact that the Nielsen-Olesen 
solution minimizes the energy 

/'OO /'OO 

E[a(r)J(r)] = £ dr = 2n 

Jo Jo 
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f + ~r(f - 1) 



dr. (6) 
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The asymptotic behaviour of the solutions for r<l (and finite A) is 



/ = /nr"- (A + 4 7 1 2 . )/ " r" +2 + ..., a = a 2 r 2 - -4tt^ +2 + •• • • (7) 
8(n + 1) 4(n + 1) 

For r » 1 we have [11] 

a-n + ^ + ..., /=( 1 + /3 S? + - (A£4) . (8) 

I 1 + ^ + - ( A > 4 > 

Equations (jj]) with boundary conditions Eq. (0) cannot be solved analytically. By 
employing a collocation method for boundary- value ordinary differential equations equipped 
with an adaptive mesh selection procedure in a compactified grid [T2], we have solved 
numerically the equations with high accuracy (global tolerance 10 -9 ) for a large range of 
values of A. In Fig. 1 we show the energy £ as a function of A for small values of A. We 
clearly see that E/n does not depend on n at A = 1, and is increasing with n for A > 1 
and decreasing with n for A < 1. That E/n does not depend on n means that the forces 
balance at A = 1, which makes it possible for solutions corresponding to vortices at arbitrary 
separation to exist [15] . 




Figure 1: Energy per vortex number E/n versus A for Nielsen-Olesen solutions with n = 
1,2,3. 

Extending the computations for larger values of A we observe numerically a logarithmic 
divergence of the energy. This is exhibited in Fig. 2. One can also see that, to leading 
order, the energy increases quadratically with the vortex number n. A detailed analysis of 
the numerical data reveals that the energy follows the following asymptotic formula 

E nUm 7T 

-^_ = _IogA + A(n) + o(l) , (9) 
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Figure 2: Energy over n 2 versus log(l + A) for Nielsen-Olesen solutions with n = 1,2,3. 
where 

( 0.47199, n = 1 
A(n) = < -2.48172, n = 2 . (10) 
{ -3.95372, n = 3 

In the next section we will prove rigorously that the energy behaves like that by performing 
an asymptotic analysis of Nielsen-Olesen solutions. 

3 Asymptotic analysis 

For our asymptotic analysis, we first split the energy Eq. ([6]) into four parts, 

r°° a' 2 

E x = 2tt / — dr , (11) 



o 



E 2 = 2ttJ -f' 2 dr , (12) 

POQ 1 

E 3 = 2tt ( n -a) 2 f 2 dr, (13) 

r°° a 

E 4 = 2ir -r(f 2 - l) 2 dr . (14) 
Jo 8 

These four contributions to the total energy correspond to the gauge field contribution (Ex), 
the Higgs dynamical contribution (E 2 and S3), and the contribution of the potential (-E4), 
respectively. 
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To study the dependence of the energy on A we differentiate Eq. ([6]) with respect to A 
and obtain 

dE_ r°° (d£_ dadS daf_d£_ d£d£_ df_d£_\ 
~dX~ J \^ + ^^ + ^^ + dXdf + ~dXdf J ) r 

= l° O ^(f 2 -l) 2 dr>0. (15) 

Here we have used integration by parts, the equations for a and / Eq. (J3]), and have assumed 
that (a' /r)(da/dX) and rf'(df/dX) vanish as r — > and as r — > oo. We see that the energy 
increases with A, and, if the energy is bounded, that / = 1 (r > 0) in the limit A — > oo. 

We will now show that the energy of the Nielsen-Olesen vortex is not bounded for 
A — > oo, but that nevertheless / will approach the singular limit / = 1 (r > 0). We start by 
considering both possibilities. If we do not have / = 1 (r > 0) in the limit, the integral in 
Eq. (fT5j) does not go to zero and E4 is at least of order A for large A. That is in contradiction 
to numerics. In Fig. 3 we exhibit E4 as a function of A. It clearly tends to n 2 7r/2 in the 
limit A —7- 00, so it is bounded in that limit. 
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Figure 3: E4 versus log(l + A) for Nielsen-Olesen solutions with n = 1, 2, 3. 

On numerical evidence, we conclude that function / tends to the singular limit / = 1 (r > 0) 
when A — > 00. (Persueing this possibility, we will later also conclude that / must tend to 
the singular limit / = 1 (r > 0) based on a series of analytic arguments alone.) The way 
this limit is approached may be understood by plotting / as a function of the scaled radial 
coordinate y/Xr. The shape of f(yXr) depends on A very slightly, reaching the profile of 
the limiting case (A = 00) very quickly, above A ~ 100. We show this fact in Fig. 4 where / 
is plotted as a function of \f~X~r for Nielsen-Olesen solutions with n = 1 and several values 
of A. The main consequence of this is that the region where / differs from 1 for large A has 
a width of order 1 / vA. 
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Figure 4: Function / as a function of y/\r for A = 0.1,1.0,10.0 and the limiting case 
A = oo. 

In the limit A — > oo, the function a satisfies the equation 

a" - -a' - a = -n (r > 0) . (16) 
r 

The general solution of this equation is 

a = n + c\rK\{r) + C2rl\{r) , (17) 

in terms of the modified Bessel functions K\ and I\. The condition for r — > oo implies 
C2 = 0, the condition a(0) = means c\ = —n, and we have 

POO 

E 3 = Tin 2 j rKlir) dr, (18) 
Jo 

which is divergent, since the integrand is of order 1/r for small r. So the energy is definitely 
not bounded in the limit A — > oo. 

That a —> n — nrKi{r) and therefore F\2 = a!{r)/r = nKo(r) as A — > oo has been shown 
before by Berger and Chen [3]. Berger and Chen study the equation for the magnetic field 
F\2- They show that the equation for F\2 linearises and is of the form 

- AF i2 (x) + F 12 (x) = 2irn5(x), (19) 

in the limit A — > oo. F\2 = nKo(r) is the solution of this equation. 

Before we derive the asymptotic behaviour of E% for large A, we calculate the A — > oo 
limit of Ei and E4.. For a = n — nrK\(r) we have 

Ex = ™ 2 [ rKl dr = ^ [r 2 (K 2 - K 2 )] ™ = ?f. (20) 




6 



Because the solution (a(r),/(r)) minimizes the energy, we have a family of functions 
(0(77-) , / (7?")) that satisfies 



-^E[a(jr),f^yr)] 



= ±- { 1 2 E 1 [a(r)]+E 2 [f(r)]+E 3 [a(r)J(r)]+7- 2 E i [f(r)]) 
7=1 a ^ 



7=1 

= 2(E 1 [a(r)]-E 4 [f(r)}) = 0, (21) 

which is a manifestation of Derrick's theorem. Therefore, both E% and £4 approach the 
finite value nn 2 /2 in the limit A — > 00, in agreement with the numerical computations (see 
Fig. 3). Using the asymptotic value of £4 in Eq. (fl~5|) . we get 

dE Tin 2 _ 7T71 2 , . . 

«=ir « E = — losA - (22) 

to leading order. Since the solution minimizes the energy, the second possibility, where / 
tends to / = 1 (r > 0), must be the one that is realized. We have already seen that the 
energy is at least of order A for large A, if / does not tend to / = 1 (r > 0). 

The logarithmic divergence of the energy of Nielsen-Olesen solutions in the limit of large 
A comes from the contribution £3, since E2 remains finite. E2 and £3 have the following 
behavior for large A: 



E 2 = n 2 A 2 (n) + o(l) , (23) 
E 3 = |n 2 logA + n 2 A 3 (n) + o(l) , (24) 

where the first three values of the functions A2(n) and A 3 (n) are 

( 0.87679, n = 1 
A 2 (n) = ! 0.32589, n = 2 , (25) 
( 0.17708, n = 3 

and 

( -3.54639, n = 1 
A 3 (n) = < -5.94920, n = 2 , (26) 
I -7.27239, n = 3 

respectively. 

Before we continue with our asymptotic analysis, we look at the variational analysis by 
Hill et al. [7] for large Ginzburg-Landau parameter. Hill et al. use the functions 

/ = l_ e -/-, a = n(l-e-^) 2 , (27) 

and minimize the energy with respect to [i and h. (From our previous discussion we know 
that /j should go to infinity and h should go to a constant as A — > 00, if there is any chance 
of approximating the correct asymptotic results.) With this ansatz the four terms of the 
energy are 

£ 1 =4 7 rn 2 / l 2 log^, E 2 = vr/4, E 3 = irn 2 G(s), E A = J^, (28) 
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where s = fj,/h and 

G{S) = lDg 2"( a + 3)8( a + l)4 • (29) 
Minimizing the energy with respect to /j, and /i leads to the equations 

h 288n 2 G"(s) h 3 G'(s) 



/i 3 89A ' n 81og(9/8)' 

For large s, G(s) = log s + log(3 4 /2 7 ) + 0(l/s 2 ), 



(30) 



M=-^+of4=V * = ; 1 +0(7). ( 31 ) 

12\/2n VVa; 2^21og(9/8) VV 



and 



^1 = —, i?2=vr/4, £ 3 = 7m 2 1 -logA + log V ^ '-L \ E 4 = —, (32) 

up to order o(l). We see that this approximation gives the correct leading terms for E\ 
and £3. Using the argument we used in Eq. (f2Tj) on the energy E(fi,h) we get E\ = E4, 
and therefore the leading term of E4 must also be correct. The 0(1) terms in E2, E% and 
the total energy E are not correct. For n = 1, e.g., the variational method gives the upper 
bound E = (vr/2) log A + 0.551, whereas the correct value is E = (vr/2) log A + 0.472, as we 
saw previously (see Eq. (|10p ). That we do not obtain the correct values is no surprise. For 
A — > 00 the function / in Eq. (|27j) goes to the step function, which is the correct asymptotic 
limit. The function a in Eq. (|27p . however, does not go to n — nrK\{r). Furthermore, 
the limit is not approached using the asymptotic expansions of solutions. The functions in 
Eq. (l27|) do not even have the correct asymptotic behaviour Eq. flHJ) for large r. 

We now give the correct asymptotic approximation for large A. Motivated by Fig. 4 and 
its interpretation, we are looking for a family of approximations with the following features: 
In the outer region, / approaches 1, and a approaches n — nrK\. In the boundary layer 
(for r < ro), / gets steeper with increasing A and the width of the layer goes to zero in the 
limit. This means that the outer approximation a = n — nrK\ extends down to r = in 
the limit A —> 00, although n — nrK\ does not have the asymptotic behavior Eq. of a, 
since 1 — rK\ = — (r/2) logr + ... for small r, i.e., the limit is singular. 

Away from the boundary layer, we look for an outer solution of the form 

/ = 1 - —f + . . . , a = n - nrKi + -a + ... (r > r ) (33) 
A A 

and find 

poo poo 

f = n 2 Kf, a = k n rK x + 2n 3 rKi / sh{s)K\{s) ds - 2n 3 rh sKf(s)ds, (34) 

J r J r 

where k n is a constant. For r > 1 the solutions are of the form Eq. ([8|) with a = w ir/2n. 
Also //A < 1 holds for r > 1/y/X. 
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In the boundary layer, a stays very small and / rises rapidly. As an approximation we 
can therefore use for the inner solution the equation 

f" + lf'- ! ^f=^f 2 - 1 )f (0<r<r Q ) (35) 

with 

/(0) = 0, /(r ) = 1 - jf(ro), (36) 

instead of using the second-order equation for / in Eq. @. The solution of this boundary 
value problem, denoted by /, has to be found numerically. Given /, we then have to solve 
the equation 

a"--a' + f 2 (n-a) = (37) 
r 

with 

a(0) = 0, a(r ) = n — nrKi(r ) + -r-a(r ). (38) 

A 

We will denote this inner function a by a. 
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Figure 5: Comparison of the numerical (exact) functions a and / for n = 1 and A = 100 
with the linear approximations a approx and /approx given by Eq. (f33|) away from the boundary 
layer and a and / in the boundary layer. 

In order to show that a good linear approximation of the functions a and / for large A 
is given by Eq. (|33p away from the boundary layer and a and / in the boundary layer, we 
compare in Fig. 5 the numerical (exact) functions a and / with the corresponding linear 
approximations a ap p r ox and / apP rox for n = 1 and A = 100. We observe that for a value of 
the location of the layer ro such that vAro ~ 20, the relative deviation of the approximation 
with respect to the exact values is of the order of 10~ 5 . This agreement improves as A is 
increased. 
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4 Conclusions 



To complete the study of the four terms which contribute to the energy, we have used 
numerical computations. The asymptotic result for the total energy Eq. (|22p . however, 
follows from a simple chain of analytic arguments, as we have seen. In contrast to Hill et 
al. [7] we make no assumptions about the class of functions to be considered. An important 
step in our chain of arguments is that in the A — > oo limit the Higgs field takes its vacuum 
value for r > 0. In this regard, the vortex behaves like the monopole [8J. The crucial 
difference is that after the Higgs field has decoupled, the energy from the interaction of the 
Higgs field and the gauge field diverges in the case of vortices, whereas it is finite in the 
case of monopoles. 
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